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We study the stability of a collisionless, relativistic, finite-strength, cylindrical layer of charged 
particles in free space by solving the linearized Vlasov-Maxwell equations and compute the power of 
the emitted electromagnetic waves. The layer is rotating in an external magnetic field parallel to the 
layer. This system is of interest to understanding the high brightness temperature of pulsars which 
cannot be explained by an incoherent radiation mechanism. Coherent synchrotron radiation has also 
been observed recently in bunch compressors used in particle accelerators. We consider equilibrium 
layers with a 'thermal' energy spread and therefore a non-zero radial thickness. The particles interact 
with their retarded electromagnetic self-fields. The effect of the betatron oscillations is retained. A 
short azimuthal wavelength instability is found which causes a modulation of the charge and current 
densities. The growth rate is found to be an increasing function of the azimuthal wavenumber, a 
decreasing function of the Lorentz factor, and proportional to the square root of the total number 
of electrons. We argue that the growth of the unstable perturbation saturates when the trapping 
frequency of electrons in the wave becomes comparable to the growth rate. Owing to this saturation 
we can predict the radiation spectrum for a given set of parameters. Our predicted brightness 
temperatures are proportional to the square of the number of particles and scale by the inverse 
five-third power of the azimuthal wavenumber which is in rough accord with the observed spectra 
of radio pulsars. 



I. INTRODUCTION 



The high brightness temperatures of the radio emis- 
sion of pulsars (Tr ^ lO^^K) implies a coherent emission 
mechanism 0,0; IS Sill and some part of the radio emis- 
sion of extragalactic jets may be coherent @. Recently, 
coherent synchrotron radiation (CSR) has been observed 
in bunch compressors j7|, |3, j9|| which are a crucial part 
of future particle accelerators. When a relativistic beam 
of electrons interacts with its own synchrotron radiation 
the beam may become modulated. If the wavelength of 
the modulation is less than the wavelength of the emit- 
ted radiation, a linear instability may occur which leads 
to exponential growth of the modulation amplitude. The 
coherent synchrotron instability of relativistic electron 
rings and beams has been investigated theoretically by 
in 111 El- Goldreich and Keeley analyzed the 
stability of a ring of monoenergetic relativistic electrons 
which were assumed to move on a circle of fixed radius. 
Electrons of the ring gain or lose energy owing to the 
tangential electromagnetic force and at the same time 
generate the electromagnetic field. 0| analyzed the sta- 
bility of a relativistic electron ring enclosed by a conduct- 
ing beam pipe in an external betatron magnetic field. A 
distribution function with a spread in the canonical mo- 
mentum was chosen for their analysis. For simplicity 
the effect of the betatron oscillations was not included in 
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their treatment. They find a resistive wall instability and 
a negative mass instability. Furthermore, they find an in- 
stability which can perturb the surface of the beam. 
analyzed the stability of a ring of relativistic electrons in 
free space including a small energy spread which gives 
a range of radii such that particles on the inner orbits 
can pass particles on outer orbits. J^J has developed a 
similar model which includes the effects of the conduct- 
ing beam pipe. Numerical simulations by fl^ show the 
burst-like nature of the coherent synchrotron radiation. 

The present work analyzes the linear stability of 
a cylindrical, collisionless, relativistic electron (or 
positron) layer or E-layer Particle densities 

in pulsar magnetospheres are very low, of order the 
Goldreich- Julian charge density ncj = ^2 • B/27rce ~ 
10"scm^3(5/;L0i2G)(i?/r)3[P(sec)]-i at radius r > R, 
where R is the stellar radius, B = 10^^i?i2 G is the sur- 
face field strength, and P is the rotational period; thus, 
the magnetospheric plasma is collisionless to an excel- 
lent approximation The particles in the layer have 
a finite 'temperature' and thus a range of radii so that 
the limitation of the Goldreich and Keeley model is over- 
come. Although we allow a spread in energies, we assume 
that it is small, so the charge layer is also thin; efficient 
radiation losses are probably sufficient to maintain rather 
low energy spreads in a pulsar magnetosphere, although 
the precise size of the spread is still not entirely certain. 
Viewed from a moving frame the E-layer is a rotating 
beam. The system is sufficiently simple that it is relevant 
to electron fiows in pulsar magnetospheres (cf. ,18]). The 
analysis involves solving the relativistic Vlasov equation 
using the full set of Maxwell's equations and computing 



the saturation amplitude due to trapping. The latter 
allows us to calculate the energy loss due to coherent 
radiation. 

In §2 we describe the considered Vlasov equilibria. The 
first type of equilibrium (a) is formed by electrons (or 
positrons) moving perpendicular to a uniform magnetic 
field in the z— direction so as to form a thin cylindrical 
layer referred to as an E-layer. The second type of equi- 
librium (b) is formed by electrons moving almost parallel 
to an external toroidal magnetic field and also forming 
a cylindrical layer. §3 describes the method of solving 
the linearized Vlasov equation which involves integrat- 
ing the perturbation force along the unperturbed orbits 
of the equilibrium. In §4, we derive the dispersion re- 
lation for linear perturbations for the case of a radially 
thin E-layer and zero wavenumber in the axial direction, 
kz = 0. We find that there is in general a short wave- 
length instability. In §5 we analyze the nonlinear satura- 
tion of the wave growth due to trapping of the electrons 
in the potential wells of the wave. This saturation al- 
lows the calculation of the actual spectrum of coherent 
synchrotron radiation. In §6, we derive the dispersion re- 
lation for linear perturbations of a thin E-layer including 
a finite axial wavenumber. The linear growth is found 
to occur only for small values of the axial wavenumber. 
The nonlinear saturation due to trapping is similar to 
that for the case where = 0. In §7 we consider the 
effect of the thickness of the layer more thoroughly and 
include the betatron oscillations. §8 discusses the appar- 
ent brightness temperatures for the saturated coherent 
synchrotron emission. §9 discusses some implications on 
particle accelerator physics. §10 gives conclusions of this 
work. 




FIG. 1: Geometry of relativistic E-layer in (a) for the case of 
a uniform external axial magnetic field, and in (b) for an ex- 
ternal toroidal magnetic field with an external radial electric 
field. 



is the electron rest mass, — e is its charge, and the units 
are such that c = 1. Here, the external magnetic field is 
assumed to be uniform, B*^ — B^z, with = r_B|/2, 
and B| > 0. Thus we have /° = f{H,P^). We consider 
the distribution function 



II. EQUILIBRIA 



f = KS{P^ - Po) exp [ - H/T] 



(3) 



A. Configuration a 

We first discuss the Vlasov equilibrium for an axisym- 
metric, long, thin cylindrical layer of relativistic electrons 
where the electron motion is almost perpendicular to 
the magnetic field. This is shown in Figure la. The 
case where the electron motion is almost parallel to the 
magnetic field is discussed below. The equilibrium has 
d/dt = 0, d/d(j) = 0, and d/dz = 0. The configura- 
tion is close to the non-neutral Astron E-layer of [l^ . 
The equilibrium distribution function /° can be taken to 
be an arbitrary non-negative function of the constants of 
motion, the Hamiltonian, 



and the canonical angular momentum, 
P4, = rp^ - erA^{r) , 



(1) 



(2) 



where = A'^ + A^ is the total (external plus self) 
vector potential, is the self electrostatic potential, rrie 



where K, Pq, and T are constants (see for example |l9j). 
The temperature T in energy units is assumed sufficiently 
small that the fractional radial thickness of the layer is 
small compared with unity. Note that a Lorentz trans- 
formation in the z— direction gives a rotating electron 
beam. 

The equations for the self-fields are 



1 d 
r dr 



dr 



= 47re / d^p f{H,P^) 



d fldirA^) 



(4) 



(5) 



where = {P^/r + eA^)/ H . 

Owing to the small radial thickness of the layer, we 
can expand radially near tq 
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where Di, D2 are the derivatives evaluated at tq, and 
Sr = r — with (Sr/ro)^ ^ 1. We choose tq so as to 
ehminate the term hnear in Sr. Thus, 



H = Ho- e$^(ro) + ^^{pl+pl + w}, 5r^) 

where ujpr is the radial betatron frequency, and 

-12^1/2 
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(7) 



(8) 



We assume 7q S> 1 and v^o > so that w^g = 1 ^ l/(27o) 
to a good approximation. The "median radius" tq is 
determined by the condition 



Di 
2Fo 



dr 
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Ho \ro / \ ' 

To a good approximation. 



dA^ 
dr 



(1 - 2C)eB| eB| 



(1 + 20 or 



(9) 



Here, 



c 



2^ 



[1 + 3C + 0(C')] 



with 



C' « 1 , 



(10) 



(11) 



is the field-reversal parameter of Christofilos. For a ra- 
dially thin E-layer of axial length L consisting of a total 
number of electrons N, the surface density of electrons is 
a = N/(2TrroL) and the surface current density is —ev^a. 
Because i3|(ro) is one-half the full change of the self- 
magnetic field across the layer, we have ( = reN/{'-fL), 
where = e^/(mc^) is the classical electron radius. No- 
tice that N, C, and 7L are invariants under a Lorentz 
transform in the z— direction. 

The radial betatron frequency ujpr is given by 



7:r2 , ,2 _ -^2 

-tin i^j3r — — 



D2 

2 



- HqC- 



dr"^ 



(12) 



Using Eq. (O gives 



1-4C > 
1-2C rl 



dr"^ 



dr"^ 



1 - 2C 



roAr7^ 



(13) 



The term oc 1/Ar is the sum of the defocusing self- 
electric force and the smaller focusing self-magnetic force. 
For the layer to be radially confined we need to have 
C < ^7172 -/^{Ar/ro). For C < 7^(Ar/ro) and < 1, 
we have copr = l/^o to a good approximation. 
The number density follows from Eq. 



n « no exp 



Sr^ 
2Ar2 



1/2 



where Ar : 



1-2C- V2AC('^o/Ar)/72 



where 



Vth 



( T 



1/2 



and 



no = 2-KKHoTrQ exp I — 



(14) 



(15) 



(16) 



As mentioned we assume the layer to be radially thin 
with (Ar/ro)^ ^ 1. Consequently, equations Q and Q 
become 



d2$« 
dr^ 

dr^ 

Thus we obtain 

c 



Aneno exp 



6r^ 
'2Ar^ 



AireriQ w^o exp 



2Ar2 



-B|(ro) AT:enoV4,oAr^/TT/2 



(17) 



(18) 



The equilibrium is thus seen to be determined by three 
parameters, 

C^ v^,^, and l/7o^ (19) 
which are all small compared with unity. 

1. Equilibrium Orbits 
From the Hamiltonian of Eq. Q we have 



d^5r 



-Luj,.Sr, ^ Sr(t')^ dr,siii[ujf3r{t'-t)+ip] , (20) 
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where r — tq = dri sin ip. For future use we express the 
orbit so that r{t' = t) = r, where (r, t) is the point of 
observation. Also, we have 



d0 _ P^ + erA4,{r) 



dt 

so that 



d4) 
dr 



Sr 



(21) 



upr or 



Sri cos[(jj^r(i' — t) + f] + Sri cos{ip) 



(22) 



where d(j)/dr\rg — —(po/rQ. For C, <C 7^(Ar/ro) and 
^ 1, we have d4>/dr\ro/uJi3r — —^/fo to a good 
approximation. Because the E-layer is uniform in the 
2— direction, 

z(t') = z+{t' - t)v^ . (23) 
The orbits are necessary for the stabihty analysis. 



B. Configuration b 

Here, we describe a Vlasov equilibrium for an axisym- 
metric, long, thin cylindrical layer of relativistic electrons 
where the electron motion is almost parallel to the mag- 
netic field. The equilibrium distribution function f'^ is 
again taken to be given by Eq. Q in terms of the 
Hamiltonian, H, and the canonical angular momentum, 
= rp^ — erA^{r), where = A'^. We make the 
same assumptions as above, 7^ ^ 1, T/{mej) <C 1, and 
Ar^/rQ <C 1. In this case there is no external field. 
Instead, we include an external toroidal magnetic field 
with corresponding vector potential A| and an ex- 
ternal electric field E*^ with potential The fields 
and E*^ correspond to the magnetic and electric fields of 
a distant, charged, current-carrying flow along the axis. 
Thus, \E^\ < \B^\. The considered external field is of 
course just one of a variety of fields which give electron 
motion almost parallel with the magnetic field. Note also 
that the distribution function is restricted in the respect 
that it docs not include a dependence on the canonical 
momentum in the z— direction Pz = rrie^Vz — eAz. 

The distribution function (3) gives Jz = so that there 
is no toroidal self magnetic field. Thus the self-potentials 
in this case are also given by equations Q and ^ . Equa- 
tions © - ® are also applicable with the replacement 
of by the total potential <I>. In place of Eq. ifTUIl we 
find 



ro 



{l-2C)eE-{ro) eE^{ro) 



(1 + 2C) 



(24) 



where ( = B^{ro)/ E^{rQ). We again have ( = reN/{-fL), 
where re = /(rac^) is the classical electron radius and 



L is the axial length of the layer. Because d'^^'^/dr'^ = 
— {l/r)d^'^/dr, the radial betatron frequency is again 
given by Eq. H13|l (with $ now the total potential) so 
that the orbits given in §2.1.1 also apply in this case. 
The electron motion is almost parallel to the magnetic 
field in that {B^/B^y ^ C{E^/B^^)^ < < 1. No- 
tice that Eq. H24|l for ro is formal in the respect that 
E^ cx 1/r. Therefore, rp is in fact arbitrary in this case. 
Because the wavelengths of the unstable modes are found 
to be small compared with tq, it may be interpreted as 
local radius of curvature of the magnetic field. 



III. LINEAR PERTURBATION 

We now consider a general perturbation of the Vlasov 
equation with /(r, p, t) = /°(r, p) -I- Sf{r, p, t). To first 
order in the perturbation amplitude Sf obeys 



d_ 

dt 



d dp d 
dr dt dp 



Sf 



D5f 
Dt 



e(5E + vX(5B) • ^ , (25) 

where (5E and (SB are the perturbations in the electric 
and magnetic fields. All scalar perturbation quantities 
are considered to have the dependencies 



F{r) exp{im(f) + ikzZ ~ icut) , 



(26) 



where the angular frequency uj is taken to have at least 
a small positive imaginary part which corresponds to a 
growing perturbation. This allows for a correct initial 
value treatment of the problem |20|. For a perturbation 
taken to vanish as t ~* —00, 

'5/(r,p,t) = 



e J' dt'|5E[r(t'),<']+v(i')x^B[r(i'),t']| • ^ , (27) 

where the integration follows the orbit [r(t'), p{t')] which 
passes through the phase-space point [r, p] at time t. For 
the considered axisymmetric equilibria, 



df^^pdf_ 

dp H dH 



rcf) 



df 



dPs 



(28) 



H 



where the partial derivatives are to be evaluated at con- 
stant and H , respectively. Thus, the right-hand side 
of Eq. H25() becomes 



IT 
dS^ 

~dr 



b S'i' - (5$) + «mv_L • SA 



df, 
dH 

df_ 

dPs 



(29) 
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where = -V5$ - dSA/dt and SB = VxJA, = 
r5A^ is the perturbation in the flux function, = 
{vr^Vz)^ and d/dt — 9/9< + v- V. We assume the Lorentz 
gauge \/-5A + d5<^/dt ^ 0. 
Evaluating Eq. (|77|l gives 



-(5$ + itj /" dt' - 5<^' + • ^A'] 



dPs 



: j dt' U'S^' - 5^' + v'l • 5 A' 



(30) 

where the prime indicates evaluation at \r(t'),t']. The 
integration is along the unperturbed particle orbit so that 
df^ /dH and df^ /dP^ are constants and can be taken 
outside the integrals. Note also that d/dt acting on a 
function of (r,t) is the same as D/Dt. 



IV. FIRST APPROXIMATION 

As a starting approximation we neglect (i) the radial 
oscillations in the orbits [(Ar/ro)^ <C 1], (ii) the self-field 
corrections to orbits proportional to C, (iii) the terms in 
6f proportional to (w^^ w (Ar/ro)^ <C 1), (iv) we 
take kz — and (v) wc assume the layer is very thin. 
Owing to approximation (iii) , we can neglect the terms cx 
vj^ • SA in Eq. (|30() in the evaluation of Sp and (5 J^. This 
is because these terms give contributions to Sf which 
are odd functions of Vr and v^. Therefore, their average 
contribution can be neglected. 

Evaluation of Eq. (|30|) gives 



dH 



m(j) 



dP^ 



H 



(31) 

where = 4>{ro). The approximations lead to a closed 
system with potentials ((5$, S"^) and sources {Sp, SJ^). 
We have 



Sp^ 



J dprdpzdP^ Sf , 



SJ^ = -e J d^P '"<t>Sf = j dprdpzdPff, v^Sf . (32) 

For the considered distribution function, Eq. 121), 
df"/dH = -f"/T. The df"/dP^ term in Eq. ETJ 
can be integrated by parts. Furthermore, note that 
dH/dP^ = (j) and d(j)/dP^ = — (0)^/-ff, which corre- 
sponds to an effective "ne gativ e mass" for the parti- 
cle's azimuthal motion [2^ l23l 01 • From the partial 
integration the small term proportional to dv^/dP^ — 
v^/{raH^) is neglected. Also note that H is not a con- 
stant when performing the integration over momenta. 



Evaluating this term by an integration by parts with a 
general function g{P^) in the integrand gives 



dP^ 



dPs 



9{P4>) 



H 



d 



K I dP^SiP^^P,)— [g{P^)e-"/^ 



K J dP^S{P^^Po)^[giP^)]e-"/^ 

^ J dP^S{P^ - P,)giP^)e-"/^^ . (33) 



That is, the integration produces an additional term 
which cancels the 1/T-term. Thus, 



dP^Sf 



where Auj — uj 
momenta gives 



dP,^- 



/° (wJ * - m(5<i>) 



(34) 



Integrating over the remaining 



{Sp, SJtjy) = (l,i'0) n{r) 



m(j>'^ {ujS'i! - mS^) 



H 



{Acoy 



(35) 



For a radially thin E-layer we may take n{r) — 
rto exp(— (5r^/2Ar^) no\/27rAr S{Sr). We comment 
on this approximation below in more detail when we 
include the radial wavenumber kr of the perturbation. 
Then equations (jA.4|l and <jA.5|) can be written as 

[<5$(ro), (5*(ro)] = [l, rov^l+AQ)] 2^^^ Z j dr Sp{r) , 

(36) 

where Z = iJm{ujrQ)Hm\^'i'f)), uj = uj/{ra(j)) and Auj = 
ALu/{m(p). Integrating Eq. I|35|) over the radial extent of 
the E-layer and cancelling out the field amplitudes gives 
the dispersion relation 



1 = 27rVo [noe^\/27rAr] Z 



m(f)^ ujroVcf,{l + Auj) — m 



H 



{Auj) 



In terms of dimensionless variables this becomes 



where Z = iJm{mCjv^)Hm\mujv^), 
Jm + iYm, and the field- reversal parameter 
Q = 47renof0Ar 1/7172/51 as given by Eq. (fTT|l . 

For m » 1, Jm{'mujv^) « (2/r7i)^/'^Ai(w) and 
Ym{mil!Vfi,) ~ —{2/m)^/^Bi{w), where Ai and Bi are 
the Airy functions and w = (771/2)2/3(^-2 _ 2Aa)) [3ll| . 

Thus we have Z = iJ^HiV « (2/m)2/3[Ai(w)Bi(w) + 

It is useful to denote Z as Zjn{w). For |wp 
1, Ai(u;) « (2V^)-iw-i/4exp(-2w3/2/3), Bi{w) w 
(0F)-iu;-i/4exp(2«;3/2/3), Z^ « (2/m)2/3/(27r|«;|i/2). 



(Ati)2 



(37) 



(38) 



r(l) 
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For |wp < 0.5, Ai{w) = ci - C2W + 0{w^) and Bi(w) = 
V^[ci + C2W + 0{w^)], where a = l/[3^/^T{2/3)] « 0.355 
and C2 = l/[3i/3p(x/3)] « 0.259. In this Hmit we have 

C2w)^]. For 



Z^iw) « (2/m)2/3[V3(c2 - cW) + 
< 1, Z,n ~ (0.347 + 0.200 i)/™^/^. 



A. Range of Validity 

We are interested in the regime where the wavelength 
of the emitted radiation is comparable to the "bunch 
length", i.e. u; « m or equivalently Auj <^ 1. How- 
ever, Eq. I2HI is only valid if Auj ^ 7^^. Since we ne- 
glected dJr and SJz we obtain from the continuity equa- 
tion SJ^ — ^^^Sp. Due to this approximation the factor 
on the right hand side can become bigger than the speed 
of light if Alo > which leads to unphysical results. 
In the latter case 6J,j, = v^6p is a better approximation. 
Fortunately, Aw <C 7^^ is the most interesting case and 
in the remainder of this paper we will always work in 
this limit. Furthermore, for the continuum approxima- 
tion to be valid the mean particle distance has to be much 
smaller than the wavelength. 
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FIG. 2: The graph shows the frequency dependence of the 
growth rate for a sample case where 7 = 30 and ^ — 0.02 
obtained from our approximations for Eq. I),'i8^ . For these 
parameters, mi « 10^ and m2 ~ 2.7 x 10^. 



B. Growth Rates 



in this regime. For case II we assume that m ^ TO2, in 
which case Eq. H39|l implies 



It will prove useful to define two characteristic val- 
ues of m: mi = (^/^^"^ and 1712 — 2^^, and therefore 
mi = (^^"^1712/2. We can obtain approximate solutions 
to Eq. (|38|) in two different cases. There may be solu- 
tions with small values of 7^ Aw, so that w ~ {m/m2)'^^^- 
In this case, Eq. (|38|l becomes a simple quadratic equa- 
tion, which can be solved for Aw. We can simplify the 
solution somewhat by changing variables to cr = 7^ Aw 
in which case Eq. (|38|) can be written in the form 



1 = -2 (^-1) 



where we have neglected a compared to one in the ap- 
proximate version of this equation. We find that 



y/-TTCZmJ^ 



(39) 



For case I let us assume that m <ti m2 , in which case Eq. 
implies 

rT~±1.121(mi/m)i/3 e'(Wi2) 
= 1.121(mi/m)i/'''(-0.2588 + 0.9659i) . (40) 

so \a\ ^ 1 for TO ^ TOi. The growth rate of the unstable 
mode is 



Wi 



(41) 



,•^1/2 3/2 

1 = ±zC'/'(m2/2TO)i/2 



,1/2 



C 



1/2^1/2, 



7 



1/2 



(42) 



note that the growth rates in cases I and II match almost 
exactly at m = m2, where \a\ ~ 

Note that TO20 is the approximate frequency of the 
peak of the single particle synchrotron radiation spec- 
trum. For more accurate results we employ a numerical 
method for solving Eq. 1381 outlined in 34] . This method 
also allows us to count the number of roots which are 
enclosed by a contour. So far we have no numerical evi- 
dence of the existence of more than one solution with a 
positive real part. The numerical results agree very well 
with our approximations even if to < toi and are shown 
in Fig. m 



C. Comparison with Goldreich and Keeley 

Goldreich and Keeley |3i] find a radiation instability in 
a thin ring of relativistic, monoenergetic, zero temper- 
ature electrons constrained to move in a circle of fixed 
radius. Under the condition 1 <C m^/"^ ^ 7 their growth 
rate is Wi « 1.160 m^^^[reN/ {'f'^ro)]^^'^ which is close to 
our growth rate with L replaced by tq. 
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V. NONLINEAR SATURATION 

Clearly the rapid exponential growth of the linear per- 
turbation can continue only for a finite time. We analyze 
this by studying the trapping of electrons in the moving 
potential wells of the perturbation. For (Ar/ro)^ ^ 1, 
the electron orbits can be treated as circular. The equa- 
tion of motion is 



dP4, 
dt 



rSF^ , SF^ = -^e[SE^ + {vxSB)^] , (43) 



where is the canonical angular momentum, where 

SF,j, = —e6E,j,Q exp{ujit) cos(m0 — Urt) , (44) 

where SE^q is the initial value of the potential, ujr = 
Re(Li;), and Ui = Im(cL;). 

For a relativistic particle in a circular orbit, 



SPcf, = me*rQS(j), where me 



r 



1 



-TOe7, (45) 



where is the "effective mass," which is nega tive, for 
the azimuthal motion of the electron ( 22, !2J| or |23l| . 
p. 68). Combining equations (|43|) and 145|) gives 



(46) 



where (p = rrup — Uft + ^tt, lot = wto exp(a;it/2), and 

ujtq = [eTO(5£'0o/(me7ro)]^/^, where ujt is termed the 
"trapping frequency." At the "bottom" of the potential 
well of the wave, sinip « ip. An electron oscillates about 
the bottom of the well with an angular frequency ~ ujt- 
This is of course a nonlinear effect of the finite wave am- 
plitude. A WKBJ solution of Eq. gives 

ip oc Wj-o^^ exp(— [.i;ii/4) sin {^{2ujto / uJi)[eyiY){uJit / 2) — 1]} . 

(47) 

The exponential growth of the linear perturbation will 
cease at the time tsat when the particle is turned around 
in the potential well. This condition corresponds to 
^T{tsat) ~ ^i- Thus, the saturation amplitude is 



\5E, 



rriej 



\^ / Wj(m)y 
\5E{tsat)\ = \6Eo\exp{uj.,tsat)- 



(48) 



where \SEsat\ 

VI. FIRST APPROXIMATION WITH fc^ / 

Here, we consider kz but keep the other approx- 
imations. Our ansatz for 5f is general enough to han- 
dle this case since it retains the biggest contribution to 
the Lorentz force in the z-direction which is of the order 
v^Br- In place of Eq. 134|l we obtain 

/ dP^Sf -e / dP^^ , X I. ^2 ' (4^) 
J J tl [io ~ rriip — KzVz) 



where we assume without loss of generality fc^ > and 
kz <C m/rQ, Lo. In place of Eq. we find 



exp(-f2/2w,\)^ ^ 



e{Lj,kz) = l+kzA{u;,kz) I 
where 

[...] . 



^27r vth 



(50) 



(lu — m(j) — kzVzY 



Here, e acts as an effective dielectric constant for the 
E-layer, and 



Aiu, kz) = TT C Z{uj, kz) 



kd 



kzT 



Z EE ^JM^' - kl)'/'] i?(p[ro(^^ - k^y^^] , (51) 

and fc0 — m/ro is the azimuthal wavenumber. The ex- 
pression for Z is from §4. An integration by parts gives 



e(M) = 1 + 



Aiu;) J 



dvz 



exp(-'i;2/2w2^) m<l)Vz/kz 



2^ < 



Vz-U 



(52) 



where the kz dependence of e and A is henceforth im- 
plicit. We can also write this equation as 



e{u) = \ + B{u) 

















Vth 


\VthJ. 



(53) 



where 



(54) 



and 



V2^ J- 



for Im(z) > 0, and 
F(z) ^ 



V2^ J- 



, exp(-xV2) rrr- 
ax h zv 27rexp 



for Im(z) < 0. The second expression for F{z) is the 
analytic continuation of the first expression to In i(z) < 
which corresponds to wave damping (see, e.g., |25| . ch. 
5). Note that terms of order A(D have been omitted. 

For m ^ 1, the factor Z = iJmiJm + iYm) can be 
expressed in terms of Airy functions in a way similar 
to that done in §5. One finds Jm[ro{Lu^ — kzY^^] « 
(2/m)i/3Ai(w), y„[ro(w2 _ ^2)1/2] ^ -(2/m)i/3Bi(u;), 



Zr~ {—] ^ Ai(u;)Bi(u;) , w ) Ai"(u;) , 
m J \ m , 



2/3 



(55) 
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tan(V|/) = k^/k^ 



FIG. 3: The figure shows the growth/damping rate uji and 
real part of the frequency Auir = ui — m<p in units of as a 
function of tan?/) = kz/k^ for m = 100 and m = 1000 for 
an E-layers with 7 = 30, C, — 0.02 and vth = 30/7'^. In the 
region of damping uJi < 0, the second expression for F[z) in 
Eq. ll^ is used. 



For tan-0 and A(2; ^ 7~^, this is the same as Eq. 
(|38|l as it should be. In general Eq. 1561 will have more 
than one unstable mode. In the remainder of this para- 
graph we will only study the largest unstable solution for 
which we recover the growth rates found in §4 in the limit 
tan?/) 0. Fig. |21 shows some sample solutions. For the 
case shown the u dependence of Z is negligible. 

General solutions of Eq. (|53|l can be obtained using the 
Newton- Raphson method ([32|, ch. 9) where an initial 
guess of {ur,Ui) gives (er,ei). This guess is incremented 
by an amount 



5Ur 




der/dur dtr/dui 


-1 








dei/dur dei/dui 







(57) 



and the process is repeated until Sr = and Si = 0. 
Fortunately, the convergence is very rapid and gives |£| < 
10^^" after a few iterations. 

Fig. 121 shows the dependence of the complex wave fre- 
quency on the tangent of the propagation angle, tamp = 
kz/kip, for a sample cases. The maximum growth rate is 
for ijj = or kz — 0- With increasing i/j the growth rate 
decreases, and for tp larger than a critical angle tpcr there 
is damping. For the damping the second expression for F 
in Eq. must be used. Roughly, we find that the crit- 
ical angle corresponds to having the wave phase velocity 
in the z— direction of the order of the thermal spread in 
this direction, that is, Ur = AtOr/kz ~ Vth- This gives 

VC f r,N 1 1 , , 
tani/'cr Y73 = 2 3r ^73 ^ — (^8) 

for mi < m < •012- Note that the dimensionless parame- 
ter which determines the cut-off at iantp cr is j'^vth- Our 
numerical calculations of ipcr give a slightly faster depen- 
dence, tan-0cr oc l/m°'^° for this range of m. Fig. 0| 
shows the m-dependence of the critical angle. It is rea- 
sonable to assume that in a particle accelerator the weak 
focusing in the z-direction sets a low limit on k^. 



where 

kz /mN2/3 / I \ 

tan ^ — — , w = [ — ^ + tan ijj — 2u tan 

k^ V 2 / \7^ / 

It is clear that e has in general a rather complicated de- 
pendence on u ^ Ur + iui and tan^. Note that the 
expression for w goes over to our earlier w for ■0 = 
noting that utanip —> Alu. 

A limit where Eq. (|53|l can be solved analytically is 
for |mP = |Aa)p/ tan^ V > Wj^, that is, for sufficiently 
small tantp. In this limit Eq. (|53ll can be expanded as 
an asymptotic series F{z) = — 1/z — — — ... 
Keeping just the first three terms of the expansion gives 



A. Nonlinear Saturation for kz ^ 

We generalize the results of §6 by including the axial 
as well as the azimuthal motion of the electrons in the 
wave. The axial equation of motion is 

me7-^ = -e [SEz + (v x SB)z] 

w —eSEzQ exp{ujit) cos(rn(/) -I- kzZ — cot) (.59) 

The approximation involves neglecting the force cx VrSB^ 
which is valid for a radially thin layer (Ar^/rg ^ 1). 
Following the development of §6, the azimuthal equation 
of motion is 

me7ro^-^ = -eSEzo cos(m(/) -I- kzZ - tot) . (60) 
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m 



FIG. 4: Critical angle for 7 = 30, C = 0.02 and Vth = 30/7^ 



Combining equations (|59|l and (|60|) gives 
e mSEth 



(1 + tan^ V) sin , 



(61) 



where f = mcj) + kzZ — ujt + |7r and twi/j = kz/k^. Be- 
cause ■0^ ^ 1 for wave growth (Eq. jSSJ), the saturation 
wave amplitude SEsat is again given by Eq. (|48|l . 



VII. THICK LAYERS INCLUDING RADIAL 
BETATRON OSCILLATIONS 

A. The Limit fc^Ar > 1 

In this section we include the small but finite radial 
thickness of the E-layer. We keep the other approxima- 



tions mentioned at the beginning of §5. In particular 
we consider fc^ = 0. In order to include the layer's ra- 
dial thickness, we consider the wave equations within the 
E-layer, 



(62) 



where 



92 



92 



is the adjoint Laplacian operator. 

Within the E-layer, we assume that the potentials can 
be written in a WKBJ expansion as 

((5$, &'^) = {K,^, if^) exp [zr7i(/) + ifcr(r — ro) — , 

(64) 

where kj. is the radial wavenumber with (fc^Ar)^ 3> 1 
{K^,K^,) are constants. This is equivalent to assuming 
that the charge density is constant between tq — Ar and 
To -l- Ar and zero elsewhere. Evaluation of the time inte- 
grals in Eq. for r = ro gives 



E 



dt'5<^' = (5$(ro,t) X 

Jn{k5ri)i^ exp{—ikij,Sri — inip) 
i(m4> + nojfjr — uj) 



(65) 



where n is an integer, k = {kf + fc^)^/^, with k^ = m/ro, 
and tantp = kr/k^. There is an analogous expression for 
the integral of S^. We have used Eq. Ij^UII for the radial 
motion with (p = assuming <C 1 and C ^ 7^(Ar/ro) 
so that LOpr — l/''o, and Eq. (|22|l for the (/)-motion with 
dipo/drlro/ujfir — — l/fQ- Using equations l|5n|l and 
the momentum space integrals (|32|l can be done to give 



n— — oo 
00 f 



H [mcj) — uj)^ 
and finally if AO <C 7^^ 



H 



(mcf) — ujy 



E 



nom^'^K^ f r'^<puj — mil — (1 — ^o)/7^ 



H 



[lo — mq 



m X ^ 



i"e-™'^-*'='**'^'J„(A:<5ri)l 

i"e-"'^-^'=*'^'''J„(fc(Sr,)' 
(to0 -t- nw^r ~ 



Fn 



(66) 



(67) 



The prime on the sums indicate that the n = term is omitted. Here, 



(68) 
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with 



X = kAr . (69) 

The 1/T terms in Eq. |S3do not cancel exactly. They 
may be neglected if 



|AcI;|2 « Fovl 



for the n — term or if 

lAwjln/m ■ 



Ad; I < 



(70) 



(71) 



for the n 7^ terms. 

For weak E-layers we have for x — > 0, Fq ^ 1 and 
Fn^o 0. In this limit we recover the results of §4. For 
X ^ 1 and 1 ^ krVQ <C fc^ro, the Gaussian factor in the 
integrand of Fn can be neglected so that one obtains 

1 2k /mT\ 

, T-— r cos — 

/2^X l^'-l ^ 2 



-inip 



even n 



Fn 



1 2ik 



2ttx \kr\ 



/"mT\ , , 

iin (^— j , odd n . (72) 



An alternative approximation for Fn can be obtained 
by using the integral representation of the Bessel func- 
tion. The remaining integral can then be computed nu- 
merically more easily. In this way we find 



Fn 



j^n^—imp 
2^ 



d6'exp [-in9^ {x^/2){k^/k- 



(73) 

For X ^ 1, 1 ^ [kfp/krY and |n| < ^/x we can approxi- 
mate sin in the exponent by a parabola at its maximum. 
We obtain 



jTL^—inip 

2=^/^r(|)Vx 



(74) 



In general F„/(i"e~*"'^) decreases as x and n increase. 
This acts to prevent the unlimited increase of the growth 
rate as m — > oo, and it ensures that the sums over n 
converge. Fig. |31shows a plot of Fq obtained by numerical 
evaluation of Eq. |731 

Within the E-layer, Eq. gives 



47re^7iom(/)^ / rg^o; — m[l — (1 — Fo)/7^ 



H 



m 



(cj — mq 

Fn 



(75) 



In terms of dimensionless variables this equation becomes 



kl = 2m^Au} 



r 



C%/2 /(l + Ac:;)(l-7-2)-[l + (Fo-l)/7^] 



(Ad})2 

Fn 



^2/1^ (w , -^n/m — A(2;)2 



(76) 




10000 



FIG. 5: Fo for vth = 0.01 and k^ro = 10* 

where k^ = i~okn k^ = r^k^p, k = tqA:, and x — kvtt- 
Notice that Eq. [Mjcan also be written as 

(5$ = C2 sin [kr{r - ro)] + C3 cos [kr{r ~ tq)] , (77) 

for ro — Ar < r < tq + Ar. For r < tq — Ar, we have 

5$ = CMior) , (78) 

since the potential must be well behaved as r — > 0. For 
r > To -f Ar, we must have 

(5$ = Ci[jm{LOr) + lYra{LOr)\ . (79) 

This combination of Bessel functions gives 5<^{r — s- oo) — > 
for the assumed conditions where Im(w) > 0. Note that 
these potentials are just the solutions of Eq. 1)62(1 in our 
approximation for 5p. The eigenvalue problem can now 
be solved by matching the boundary conditions. How- 
ever, we have not solved the full eigenvalue problem. In- 
stead we consider unstable solutions with the restriction 
that krAr ^ 1. Under this condition we can interpret 
Eq. [THI as a local dispersion relation. Unstable modes 
found from Eq. 1761 will need a slight correction in order 
to satisfy the boundary conditions. 

We expect that Eq. [THIhas solutions near each betatron 
resonance at AQ = ±n/m. This is a familiar concept in 
the treatment of resonances in storage rings (cf. [23| or 
[23 ). We extract each solution by summing over a single 
value of n and —n only and obtain from Eq. 1761 for the 
case n 7^ and Aw ^ 7"^ 



_ _ l^VthV'^ {kl + 771^7 2) 



F^r 



Fn 



i^ + ACj) 



i^-A^S" 



(80) 



Thus, 



Ad, « ^" ± - 

^ m 
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for sufficiently big S, i.e. we expect the imaginary part 
of A(D to be negligible for the n ^ modes. Despite a lot 
of effort we were not able to prove this statement under 
more relaxed conditions. 

We can easily find an analytic solution of Eq. [TSlfor the 
case where the n = term is dominant. If |Aa)| <C 1/7^ 
and |A(D| ^ Fof-f'^, we obtain 



Aw = ±- 



(82) 



The dependence of the growth rate on fc,. becomes sig- 
nificant when ^"^k^/rn? is comparable to unity. For 
TO ^ TO2 ^ 7'^, we see that this happens when 
~ 1' which involves the combination 

"i'^vth again. 

The growth rate of Eq. |S21is proportional to y^. This 
implies from §5 that the emitted power scales as the 
square of the number of particles in the E-layer which 
corresponds to coherent radiation. Sample results are 
shown in Fig. El We conclude that the main effect of the 
betatron oscillations is an indirect one. The radial mo- 
tion itself is unimportant for the interaction. However, 
the influence of the radial motion on the time dependence 
of the azimuthal angle </) of a particle is important since 
a shift in </) can take the particle out of coherence with 
the wave. This effect is accounted for by Fq. 



B. Qualitative Analysis of the Effect of the 
Betatron Motion 

Let us suppose that vth ^ 1/7^: and that |A(I;| is not 
necessarily small compared with Vth (We can still as- 
sume |Aci)| ^ 1 without requiring the more restrictive 
condition 7^|A(i| <C 1.). The key effect of the beta- 
tron oscillations is to "wash out" the phase coherence 
of the response within the layer; for a cold layer, all 
orbiting particles move in "lock step" , which is partic- 
ularly favorable for a bunching instability. Let us sup- 
pose that |A(i| has a real part that is substantially larger 
than 1/7^. The response in the layer scales as an Airy 
function with argument w{\ +C) where |^| < Vth- The 
phase accumulated across the layer thickness ~ VthTo is 



3/2 



if A(Dr < Vth and 77 ~ mv^^^ {AQr / vth)^^'^ if 
AQr 3> Vth- Large 77 ought to imply substantial deco- 
herence of the response in the layer. We see that this is 
likely irrespective of the value of Auir/vth provided that 
TO :$> v~i^^'^, i.e. for TO/7^ ^ (7^wth)~^^^- At large val- 
ues of j'^Vth, phase smearing should suffice to suppress - 
if not eliminate - the bunching instability at frequencies 
near the synchrotron peak. Moreover, if j'^Vth ^ C~^j 
the instability should be suppressed over the entire range 
TO > (^^/^7'^ for which we found unstable modes in §4. 
Large Aujr/vth would merely accentuate the smearing. 
At a given value of to, we see that Aujr > {^"^^th)'^^ 
i.e. 7^Aci;r ^ {'m/^^)~^{'-^'^vth)~'^ suffices for large phase 
decoherence in the layer. 




m 



FIG. 6; Growth rates in the limit krVth ^ 1 for our reference 
case 7 — 30, C, — 0.02 and Vth = 1/7^ ^-^d various values of 
kr- The line proportional to m^^''* is shown for comparison. 



C. The Limit fc^ Ar < 1 

In order to determine the lowest allowed value for 
and the highest possible growth rate the full eigenvalue 
problem has to be solved. We estimate the result by eval- 
uating Eq. I A. 41 in the thin approximation again. Looking 
at Eq. I A. 41 and replacing the Bessel functions by their 
Airy function approximations for the case m 3> toi and 
TO <C TO2 we see that the thin approximation is justified if 
krVth *C 1 and rn^/^Vth <C 1. It starts to fail completely if 
mF'^'^vth ^ 1, i.e. once we start integrating over the oscil- 
lating and/or the exponentially damped/increasing part 
of the Airy function, which implies we would like to have 
TO^/^jAwl <C \/Fo with |Ad;p <C Fovff^ from the previous 
paragraph. However, for real values of kr we expect that 
the thin approximation will still give us an upper bound 
of the growth rate because it is easier to maintain coher- 
ence if all the radiation is emitted from the same orbit. 
With Eq. EZIwe obtain in the limit Alo ^ 7"^ 



-n C Z 



■Fo7 ^ 
(Ac;;)2 



(83) 



The growth rates can be found as before. For m 3> toi 
we obtain 



1/2^2/3, 



and 



1.083C'/^TO 

7 



>l/2^1/2. 



(84) 



^» - 172 V ^0 



7 

for m ^ TO2, i.e. there is an additional factor of ^/Fq- 
The results for our reference case are plotted in Fig. [7| 
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m 



FIG. 7: Solutions of the dispersion relation in the presence 
of betatron oscillations in the limit krVth ^1,7 = 30, ( — 
0.02. Points which do not satisfy the inequalities m 3> 1, 



m^^^vth < 1 and 



Aojp < Fov^f^ are plotted in gray. 



the complicated nature of the dispersion relation we have 
not done this yet. Note that the thin approximation will 
suppress certain modes, e.g. the negative mass instability 
cannot be expected to be present with the fields having 
been evaluated at one radius only, cf. p^. 



VIII. SPECTRUM OF COHERENT RADIATION 

Having computed the growth rate and the saturation 
amplitude, the radiated power can now be calculated. 
Starting from Eq. IjA.lOjl we now have 



Pm = ■T:L^rr, \5J, 



00 



(85) 



where ^ = r/ro and the integration is over the thickness 
of the layer. The Bessel function can be expressed ap- 
proximately in term of an Airy function as done before. 
We take the linear approximation to the Airy function as 
discussed previously, and this gives 



ttLuj 



.4 „2 



4/3 

— I |(5J0o| 



(86) 




where C2 ~ 0.259. This is valid for sufficiently big values 
of 7 and low m. The largest values occur for krVth ^ 
1, where this quantity is simply 4t;^^. This is enough 
motivation for us to work in this limit. Thus, 



/ o \ 4/3 

y2V^,\SJ,of [^-^) ■ (87) 



Pm < 2'KLLUr', 



Because we calculated our growth rates in the thin ap- 
proximation for A: « fc^ it is consistent to use Scj) = 
4:T:'^vthV^^ ZroSJ^Q. Furthermore, we set uj mcf). This 
is consistent even for large growth rates since the expo- 
nential growth has stopped. With our expression for the 
saturation amplitude we obtain 



P < 



Lclvlml 76 
87r3roe2 |Z|2 



4/3 



1 



uji(m) 



FIG. 8: Fo as a function of m for 7 = 30, C, = 0.02 and various 
vth and the squared ratio of the growth rates from Fig. Inland 



Fig. (dash-triple-dotted) 



which were computed numerically. In Fig. |Slthe function 
Fq is plotted which we compare with the squared ratio 
of our new growth rates to the ones evaluated previously 
without betatron oscillations. 

We could also study the effect of the non-zero thickness 
alone without betatron oscillations setting Fq ^ 1 and 
Fn-tQ = and solving the full eigenvalue problem. Due to 



Since the number of particles N is proportional to C and 
the growth rates are proportional to for > "^1 
the radiated power scales like N"^ . This suggests that 
the emitted radiation is coherent. In Fig. |51we plotted 
the radiated power in arbitrary units having evaluated 
Fq numerically. For large m the curve scales as m~^/^. 
Analytically we obtain with our second approximation 
for Fo the scaling m~'^(m2/3/m^/^)^ — m^^^^. With 
« 4ct we obtain 



P^<3.71xl0iVm-^-P^ ^ 



(89) 
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m 

FIG. 9: Radiated power {m-'^{uji/<j))^) for 7 = 30 and C = 
0.02 in arbitrary units. The straight line is proportional to 
m"^''^ and is shown for comparison. Points which do not 
satisfy the inequalities m ^ 1, m^^^vth < 1 and |A£l}|^ < 
FoVth are plotted in gray. 



IX. APPLICATIONS IN ACCELERATOR 
PHYSICS 

The next-generation linear collider requires a beam 
with very short bunches and low emittance. That is, the 
beam must occupy a very small volume in phase space. 
The emittance of the pre-accelerated beam is reduced in 
a damping ring which is operated with longer bunches to 
avoid certain instabilities. The bunch length has to be de- 
creased in a so-called bunch compressor before the beam 
can be injected into the linear collider. A bunch compres- 
sor consists of an accelerating part and an arc section. 
Since the bunch lengths of the proposed linear colliders 
are in the order of the wavelength of the synchrotron 
radiation which is being radiated in the arc section, in- 
stabilities due to coherent synchrotron have to be taken 
seriously. For a design energy of 2 GeV and 7 x 10^^ elec- 
trons per 100 fim our dimensionless quantities become 
7 = 4000 and C, — 0.08 27]. Our qualitative analysis of 
the betatron motion suggests that CSR is suppressed for 
a minimum energy spread of vth > C~^7~^ = 12.67"^. 



X. DISCUSSION AND CONCLUSIONS 



A. Brightness Temperatures 

We consider the brightness temperatures Tb for condi- 
tions relevant to the radio emissions of pulsars. Using the 
Rayleigh- Jeans formula B^, = 2kBTB{v / cY for the radi- 
ated power per unit area per sterradian at a frequency 
V = m(j)/2TT gives 

2kBTB{iy/cfAAn = 27rP,„/0 

Tb < 4.5 X 10^1 - . L7 (90) 
m V / 

where fc^ is Boltzmann's constant and A — ^-kt^L is the 
area of the E-layer. The solid angle of the source seen 
by a distant observer has been computed in appendix 
A and its value is Afi = 47r^ro/(mL). It is assumed 
that the angular size of the source is small such that 
that radiation from the top and the bottom emitted at 
an angle d with respect to the normal is received by the 
observer at the same position. For the sample values 7 = 
1000, C = 0.08, Vth = 0.047-2, L = 100 km and m = mi 
our model predicts a maximum brightness temperature of 
Tb w 2 X lO^^K. According to our results from previous 
sections there may be degeneracy from modes with non- 
zero axial wavenumbers kz. It is reasonable to assume 
that this will increase the brightness temperature by a 
factor in the order of m tamper- Beaming along the z- 
axis may increase the brightness temperature and the 
observed frequency even further. 



This work has studied the stability of a collision- 
less, relativistic, finite-strength, cylindrical electron (or 
positron) layer by solving the Vlasov and Maxwell equa- 
tions. This system is of interest to understanding the 
high brightness temperature coherent synchrotron radio 
emission of pulsars and the coherent synchrotron radi- 
ation observed in particle accelerators. The considered 
equilibrium layers have a finite 'temperature' and there- 
fore a finite radial thickness. The electrons are consid- 
ered to move either almost perpendicular to a uniform 
external magnetic field or almost parallel to an external 
toroidal magnetic field. A short wavelength instability 
is found which causes an exponential growth an initial 
perturbation of the charge and current densities. The 
periodicity of these enhancements can lead to coherent 
emission of synchrotron radiation. Neglecting betatron 
oscillations we obtain an expression for the growth rate 
which is similar to the one found by Goldreich and Kee- 
ley |3] if the thermal energy spread is sufficiently small. 
The growth rate increases monotonically approximately 
as where m is the azimuthal mode number which 

is proportional to the frequency of the radiation. With 
the radial betatron oscillations included, the growth rate 
varies as m^^^ over a significant range before it begins to 
decrease. 

We argue that the growth of the unstable perturbation 
saturates when the trapping frequency of electrons in the 
wave becomes comparable to the growth rate. Owing to 
this saturation we can predict the radiation spectrum 
for a given set of parameters. For the realistic case in- 
cluding radial betatron oscillations we find a radiation 
spectrum proportional to m~^^^. This result is in rough 
agreement with observations of radio pulsars The 
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power is also proportional to the square of the number of 
particles which indicates that the radiation is coherent. 
Numerical simulations of electron rings based on the fully 
relativistic, electromagnetic particle-in-cell code OOPIC 
[sof recovers the main scalings found here. 
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APPENDIX: GREEN'S FUNCTION 

The Green's function for the potentials give 
5^{r, t)^ j dt'd^r' G(r -r',t~ t') 5p{r' , t') . 
5A{t, t)= j dt'd^r' G(r - r', t - t') 53 [r' , t') , (A.l) 
where 

- ^) G{v,t) = -47r5(t)<5(r) , G(k,c.) = 

^, N 47r /■ , /" exp(ik • r — zwi) 
G(r,t) -J^.jj^j , (A.2) 

where G is the Fourier transform of the Green's function. 
The "C" on the integral indicates an cj— integration par- 
allel to but above the real axis, lm{u}) > 0, so as to give 
the retarded Green's function. 

Because of the assumed dependences of Eq. we 
have for the electric potential, 

/"OO pOC /'27T 

S^ujnikAr)^'^ / r'dr' / KdK / da (5pa;mfc, (j'') [■•] 
Jo JQ JQ 
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. f Ijlf J Jm{Kr)Jm{Kr') , 

= An r dr ndn. — — — dp^rnk, [r ) , 

Jo Jo 1^ - - ki) 

(A.3) 

where 

r 1 _ exp{ima)Jo{K[r^ + (r')"^ — 2rr' cosa]^^^} 

where = + ky. Because to has a positive imaginary 
part, this solution corresponds to the retarded field. Also 
because Im(cij) > 0, the k— integration can be done by a 
contour integration as discussed in [s^l which gives 

poo 

5*^™feJ0 = 27r2i / r'dr'J„,{kr<)Hil\kry)6p^^kAr'), 
Jo 

(A.4) 

where k = {lj^ — k^y^^, where r< (r>) is the lesser 

(greater) of (r, r'), and where Hm\x) — Jm{x) + iYm{x) 
is the Hankel function of the first kind. From the Lorcntz 
gauge condition 

j^^mfe.(^) = rM^"''^ = (1 + Ac;)) (5$""'^^(r) 

(A.5) 

Equations (|A.4ll and ljA.5|) are useful in subsequent cal- 
culations. 

To determine the total synchrotron radiation from the 
E-layer it is sufficient to calculate 5 A at a large distance 
from the E-layer. We assume that the E-layer has a finite 
axial length and exists between — L/2 < z < L/2. Thus 
we evaluate SA in a spherical coordinate system R = 
{R, 9, (p) at a distance L. The retarded solution is 

(5A(R) = -^ f d^r' S3(r',t-\Il-r'\ 



exp(za;i?) 
R 



5J(r')exp 



im(j)' + ikzz' — icL)(f + R • r') 



(A.6) 

(see, e.g. ch. 9 of 21]). The source point is at 
{x' = r' coscj)', y' — r' s'mcj)', z'). The observation point 
is taken to be at (a; = 0, y = RsinO, z — Rcosd). 
Consequently, R • r' = r' sin 9 sin (/>' + z' cos 9. The phase 
factor exp(iu;i?) does not affect the radiated power and 
is henceforth dropped. 

For the cases where 5 is the dominant component of 
the current-density perturbation we have 



R 



r'dr'dct)' 



m,4>. Carrying out the cf)' integration in equation ljA.7|) 
gives 



S{0) 
R 



r' dr' SJ^ir') [. 



(A.9) 



where 



iJln{ujr' sin 9), 



cor' sin 9 



J„i {ujr' sin 9) 



and where the prime on the Bessel function indicates its 
derivative with respect to its argument. The radiated 
power per unit solid angle is 



dn 



i?2 



R^ 
Stt 



R^ 

Stt 



Ikx (5A"P 



{\SA-\'+cosH\SA-\') , (A.IO) 



where k = wR is the far field wavevector. 

For a radially thin E-layer, (Ar/ro)^ ^ 1, equations 
(|A.9p and (|A.10|I give 



dn ~ Stt 
SH9) 



r'dr' 6J^{r')u)J^{ujrosin9) 



+ 



Stt 



r'dr' SJ^{r') Jm(ti;ro sin0) 



ro 



(A.ll) 



The factor within the curly brackets is the same as that 
for the radiation pattern of a single charged particle (see 
ch. 9 of 0). 

The factor S'^{9) in Eq. (jA.lip tightly constrains the 
radiation to be in the direction 0* — cos^^ {kz/tu) if 
the angular width of S'^{9), the half-power half- width 
A01/2 ~ T^/i'^L), is small compared with the angular 
spread of the single particle synchrotron radiation, I/7, 
which is the angular width due to the Bessel function 
terms in Eq. IjA.lip . This corresponds to E-layers with 
L 3> n^juj — nrQ^/m. For L ^ rg, we need m 3> ttj, 
which is satisfied by the spectra discussed later in §8. In 
this case, Eq. (|A.11|) can be integrated over the solid 
angle to give 



Pu. = 



ttL sin ( 



2uj 



r'dr' (5j0(r')ijj J^j(cjro sin6'*) 



r dr dJ^[r ) (wro sm t/* j 



ro 



(A.12) 



5J^{r') exp{im(f>' — iujr' sin 9 sin (/)') , 



(A.7) 



One limit of interest of Eq. IA.12I is that where kz = 
so that 0* — tt/2 and 



where 



S{9) = L ^^^[(^^ - wcos6i)L/2] 



(A., 



{kz - LOCOS 9)L/2 
is a structure function accounting for the finite axial 
length of the E-layer, and lo superscript indicates to = 



P„ 



nrnv^L 
2ro 



r'dr' 6J^{r')J'^{uoro) 



, (A.13) 



where we have set lo — > mcf). The total radiated power is 



